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Abstract 

The dynamics of models described by a one-dimensional discrete nonlinear 
Schrodinger equation is studied. The nonlinearity in these models appears due 



> 

I/"} ■ to the coupling of the electronic motion to optical oscillators which are treated 

in ; 

in adiabatic approximation. First, various sizes of nonlinear cluster embedded 
\q • in an infinite linear chain are considered. The initial excitation is applied 



either at the end-site or at the middle-site of the cluster. In both the cases we 
obtain two kinds of transition: (i) a cluster-trapping transition and (ii) a self- 
trapping transition. The dynamics of the quasiparticle with the end-site initial 
O ■ excitation are found to exhibit, (i) a sharp self-trapping transition, (ii) an 



amplitude-transition in the site-probabilities and (iii) propagating soliton-like 



waves in large clusters. Ballistic propagation is observed in random nonlinear 
systems. The effect of nonlinear impurities on the super diffusive behavior of 
random-dimer model is also studied. 

PACS numbers : 52.35.Nx, 63.20.Ls, 61.43-j 
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I. INTRODUCTION 



Strong interaction with the lattice vibrations is one of the basic mechanisms influencing 
the transport of quasiparticles such as electrons or exitons in solids. The consequences have 
been investigated employing different methods. The recent approach to this problem is 
based on nonlinear equations. 0-0] One of the simple models with varieties of applications 
in different areas is one-dimensional discrete nonlinear Schrodinger equation: 0-1^ 

dc m 2 

^ ^ ^(Cm+l C-rn-- .1 ) (^m Xm|Cm| ) C-m • ) 

Here c m (t) is the probability amplitude of the quasiparticle at site m at time t, V is the 
nearest-neighbor transfer matrix element, e m and Xm are on-site energy and nonlinearity 
strength of the m-th site respectively. Without any loss of generality we assume V = 1. 
The Eq. ([I]) arises in the general problem of polaron formation due to the coupling of 
quasiparticles with optical oscillators in adiabatic approximation. The simple form of the 
Eq. ([]]) with e m = and Xm = X f° r an m nas been studied numerically since long back. 
|| However, for a two-site system which is called the nonlinear adiabatic quantum dimer, 
the self-trapping transition [|]|| occurs at a critical value of nonlinearity arbitrary initial 
conditions. [|6|-|9|,|TT|1 The applications of the nonlinear dimer analysis have been made to 
several experimental situations. They are, neutron scattering off hydrogen atoms trapped 
at the impurity sites in metals, [|7| fluorescence depolarization, || muon spin relaxation, 
I"T| nonlinear optical response of superlattices [|l3l etc. The self-trapping transition also 



occurs in the extended nonlinear systems. [12| A possible application is the trapping of 
hydrogen ions around the oxygen atoms in metal hydrides. [I4j All these studies have been 
performed for finite number of nonlinear sites by assuming that the quasiparticle is localized 
within the nonlinear sites. However, the effect of nonlinear sites embedded in a host lattice 
on the dynamics of quasiparticles has been hardly studied in spite of its importance in 
real systems. Dunlap et al [|1^] studied the self-trapping transition at a single nonlinear 
impurity embedded in a host lattice. Chen et al |11 studied the time averaged probability 
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at the initial occupation site in an infinite linear chain containing one or many nonlinear 
impurities. In this study the adiabatic assumption has been removed. They have also studied 
the adiabatic case albeit not in details. So, in this paper we plan to study first the dynamics 
of a quasiparticle in an infinite linear chain containing adiabatic Holstein type impurities. 
|]J We use two different kinds of initial conditions. The initial excitation is applied either 
at the end-site or at the middle-site of the cluster of the impurities. 

If we consider randomness in site energies e m and Xm = for all m in Eq. (||) Anderson 
theory [[17] predicts that the particle will remain localized within a finite region of the chain 
after a sufficient time. So, it is important to investigate the dynamics of (i) random non- 
linear systems (randomness in the nonlinearity parameter) and (ii) systems where disorder 
in site energies and nonlinearity coexist. Regarding the first question Molina and Tsironis 
18] showed the ballistic propagation of the untrapped electronic fraction in nonlinear ran- 



dom binary alloy. In this paper we study the transport properties of completely random 



nonlinear systems. Feddersen 19 1 has studied the effect of nonlinearity on the Anderson lo- 



calization. Shepelyansky |2DJ has obtained subdiffusive behavior in on-site energy disordered 
systems only when the nonlinearity parameter exceeds a critical value. It is well known that 
super diffusive behavior is obtained in the random-dimer model (RDM). |HJ The RDM is 
characterized by a set of nonscattered states around the dimer energy. This leads to the 
super diffusive behavior of the mean square displacement of a particle. We study here the 
effect of nonlinearity on the superdiffusive behavior of the RDM. 

The organization of the paper is as follows. In Sec. II, we study the dynamics of the 
quasiparticle in a cluster of nonlinear sites embedded in a lattice. The initial excitation is 
applied either at the end-site or at the middle-site of the cluster. In Sec. Ill, we study the 
dynamics of different kinds of random systems. We end this article by summarizing our 
main results. 

II. CLUSTER OF NONLINEAR IMPURITIES EMBEDDED IN A LATTICE 
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A. Initial excitation at the end of the cluster 



We consider a system containing a cluster of n number of nonlinear impurity sites of 
equal strength x embedded in a host lattice. All the site energies are assumed to be zero. 
The initial excitation is applied at left end-site of the nonlinear cluster. We call this zeroth 
site. The sites on the left and the right of the initial occupation site are numbered as 
m = —1, —2, —3, • • • and m — 1, 2, 3, • • • respectively. We first study here the time averaged 
probability of the nonlinear sites. For m-th site it is defined as 



with |c m (0)| 2 = 5 m fl. Here, |c m (t)| 2 is the probability of the quasiparticle at m-th site at 
time t. We solve the first order coupled nonlinear differential equations numerically by using 
4-th order Runge-Kutte method. The system is taken as self-expanding lattice to avoid the 
boundary effect. For time averaging we have taken T = 200 with interval AT = 0.01. The 
accuracy of the numerical integration is checked through total probability. Here, we consider 
the cases for n = 2,3,4,5 and asymptotically large value of n {n — > oo). For n — > oo we 
mean that the system contains two semi-infinite chains. The perfect chain without any 
nonlinearity is connected to the other one which is a perfect nonlinear chain. The initial 
excitation is applied at the junction where the nonlinear impurity exists. In Fig. 1 we have 
plotted the time averaged probability at the initial occupation site, < Po > as a function of x 
for different values of n. In all these cases we find that < P > starts increasing significantly 
from x — Xcri- F° r n — 2 the value of Xcri is ~ 2.8. In general the value of Xcri increases 
with increasing the size of the nonlinear cluster. A sharp transition in < P > occurs at 
Xcr2 ~ 4.23 for all values of n. In the region between Xcri and Xcr2 we obtain fluctuations 
in < P >. For a better understanding of this behavior we study next the time averaged 
probability of the unoccupied nonlinear sites. 

In Fig. 2 we plotted the time averaged probability of the initially unoccupied nonlinear 
site (< Pi >) of the dimer as a function of x- For comparison we have also plotted < P >. 




(2) 
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When the nonlinearity strength \ exceeds Xcri we observe both < P > and < Pi > increase 
with increasing x an d their values are almost equal. This implies that the particle oscillates 
with a finite probability among the dimer sites. This partial localization or trapping within 
the cluster can be understood from the following way. At t = the energy level of the sites 
m = —1, and 1 are 0, — x and respectively. With increasing the time the site-probability 
of m = decreases and of m — 1 increases. Consequently, the energy level of the sites m = 
and 1 moves upward and downward from the original position respectively. Thus, the energy 
gap between the sites m = and 1 becomes smaller than the gap between m — — 1 and 0. 
So, the initially localized particle at site m = favors the nearest-neighbor nonlinear site 
(i.e. m — 1) and the energy level of that site decreases. At the same time the energy gap 
between the sites m — 1 and 2 increases. Thus the particle feels a quantum well and it 
oscillates within the well with a finite probability. Of course, some probability will escape 
along the leads in both directions of the dimer. Now the leakage of the probability through 
the leads reduces with the increase of x due to the increase in the energy gaps between 
the sites m = and —1 and between m — 1 and 2. When x attains a critical value, say 
Xcri, these gaps become sufficiently large to trap the particle within the cluster. So, Xcri 
marks the onset of the cluster-trapping transition of the particle. With a further increase 
in X-, the competition among the nonlinear sites to trap the particle starts. Depending on 
the strength of the nonlinearity the particle is preferentially trapped either at the initially 
occupied site or at the unoccupied site (see the inset of Fig. 2). When x is just below Xcr2, 
the value of < Pi > is much larger than < P >. But, when x crosses Xcr2 we find sharp 
fall in < Pi > and an increase in < P >. For further increase of Xi < Po > increases and 
< Pi > decreases gradually. We do not obtain any further transition. So, Xcr2 is called the 
critical value of chi for self-trapping transition. 

To understand the behavior in the fluctuation regime of time averaged site-probabilities 
we study their temporal behavior. This is shown in Fig. 3. We find both the site-probabilities 
oscillate against each other initially for some time. For a given value of x-, the time period of 
the oscillation increases and the amplitude decreases with increasing time. But after a few 
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oscillations, we observe a transition where the amplitude of the two oscillations decreases sud- 
denly and the phase is just opposite to each other (see Fig. 3(a)). The amplitude-transition 
in the site-probabilities always occurs simultaneously. Furthermore, the transition occurs 
from the peak of the oscillation in one case and from dip in the other one. Consequently, 
one of the site-probabilities oscillates with a more mean probability than the other. Thence- 
forth, the amplitude of the oscillation of the probabilities decreases with time. This kind 
of transition is obtained in the nonadiabatic nonlinear quantum dimer problem in the pres- 



ence of rapid vibrational relaxation caused by the damping in the lattice vibration. [|22|] Two 
transitions are observed, static transition at x — 2V^ || and a dynamic transition at \ = 4V \ 
H The static transition is governed by the relaxation term of the lattice vibration. Here, it 
seems that the leads at both the ends of the nonlinear dimer introduce effectively a damping 
term in the lattice vibration. This effect appears through the escape probability from the 
dimer cluster. However, the main difference here is that the transition occurs at different 
values of \. The transition time as well as the number of the dynamical adiabatic dimer type 
oscillations decreases with increasing x (compare Fig. 3(a) and 3(b)). Furthermore, near 
Xcr2 the number of the dynamical adiabatic dimer kind oscillations does not reduce for a wide 
range of \. Consequently, we do not find any fluctuation in time averaged probability in this 
region (see the region 3.87 < x < 4.23 of the inset of Fig. 2). The amplitude-transition in 
the site-probabilities occurs after half of the period of the oscillation (see Fig. 3(c)). Thus, 
just below Xcr2, < Pi > is found to be much larger than < P >. When x just crosses Xcr2 
the transition occurs at time which is even less than the half period (see Fig. 3(d)). As there 
is no dynamical adiabatic dimer type oscillation we do not find any further transition in the 



amplitude of the site-probabilities with the increase of x- So, in contrast to the Ref. fL6| , we 
find that the quasiparticle recognizes both the nonlinear impurities just above Xcri and for 
further increase of Xi amplitude-transition of the site-probabilities occurs. Furthermore, we 
obtain damped oscillation in the site-probabilities for x > Xcr2- This has to be contrasted 
with the regular oscillation above x = 4 in an isolated nonlinear dimer. |]] We also study 
the time averaged probability at the nonlinear sites for n = 3, 4 and 5. In all these cases we 
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observe the same behavior as observed in the dimer embedded in a host lattice. 

For a relatively large size nonlinear cluster (e.g. n = 30) the partial localization of the 
particle is found to occur at different region of the cluster. We study the particle propagation 
in a lattice for the case of n = 30 for different values of x- For small value of x (x < 3.1) 
we obtain the derealization behavior of the particle. For further increase of x we find a 
soliton-like wave, extended over a few sites, moves through the nonlinear cluster. Thus, we 
obtain an oscillation of the soliton-like wave in the cluster. This is shown for x = 3.51 in Fig. 
4. It should be noted that the maximum probability is found within this soliton-like wave. 
The time period of the oscillation of the wave within the cluster is increased with time. This 
indicates that if we increase the time, the oscillation will cease and the soliton-like wave will 
be localized within a few sites of the nonlinear cluster. This is exactly obtained in Fig. 4. 
The number of oscillation decreases before the localization with increasing x- The position 
of the localization of the wave, however, depends on the value of x but the pattern is not 
discerned here. For further increase of x we do not find any oscillation. It moves along 
the cluster and gets trapped within a few sites of nonlinear cluster. With increasing the 
value of x the trapping region in general moves towards the initial excitation site and we 
obtain a sharp self-trapping transition at Xcr2 ~ 4.23. It should be noted that the width 
and peak-value of the soliton-like wave decreases and increases respectively with increasing 
the value of x- Though we are not able to probe the value of x where the soliton-like wave 
starts to form but it is within x — 3.1 and 3.2. It should also be noted that the movement 
of the soliton-like wave is obtained only when the size of the nonlinear cluster is much larger 
than the width of the wave. 

In asymptotic limit (i.e. n — > oo) we also obtain the formation of soliton-like wave 
extended over a few sites in the nonlinear cluster in the system. Here, we do not find any 
oscillation of the wave as obtained in the case of n = 30. For lower values of x the wave 
moves along the lattice but finally it is trapped in a region as shown in Fig. 5. If we increase 
the value of x the localization regime of the wave moves towards the initially occupied site 
and thus obtain a sharp self-trapping transition at Xcr2 = 4.23. The width and the peak- 
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value of the wave decreases and increases respectively with increasing x- However, we did 
not probe the region of \ where the soliton-like wave starts to form. 



B. Initial excitation at the middle of the cluster 

We study here the same system but the initial excitation is given at the middle site 
(m = 0) of the cluster. The cluster contains odd number of sites. As the system is symmetric 
around the initial occupation site we do not find the asymmetric probability distribution. So, 
the properties in this system should be different from the earlier cases. The time averaged 
probability at the initial excitation site is shown in Fig. 6 for different values x an d n. In 
case of single nonlinear impurity system we obtain the self-trapping transition at Xcri ~ 3.2. 



15| , |T6H For higher values of n we find < Pq > increases significantly from x — Xcri and it 
characterizes the cluster-trapping transition. The value of Xcri for n = 3 is ~ 2.4 which is 
much less than the self-trapping transition value of x f° r n — 1. The value of Xcri hi general 
increases with increasing the size of cluster. In asymptotic limit (i.e. n — > oo) the value of 
the transition point is Xcri ~ 3.5. |23[ We further study the time averaged probability of the 



neighboring nonlinear sites of the zeroth site for n = 5 (see Fig. 7). We find that the time 
averaged probability of the other nonlinear sites also starts increasing from Xcri- The value 
of < P m > decreases as we go away from the initial occupation site i.e. as \m\ increases. 
This means that beyond Xcri the particle lies within a few sites of the cluster with center at 
the initial excitation site. Again, as both sides of the zeroth site contains nonlinearity the 
particle is attracted by the nonlinear sites in both directions. Thus, with increasing the size 
of the cluster we find that < P > decreases and consequently the value of Xcri increases. 
Beyond Xcri we find that the time averaged probability of the nonlinear sites (except the 
zeroth site) first increases and then decreases with increasing the value of x- But < Pq > 
gradually increases with increasing the value of x- Thus, we obtain the localization of the 
particle at m = with maximum probability which is called self-trapping. However, in this 
case we do not obtain any sharp self-trapping transition as seen in the previous case. In 
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the study of particle propagation we observe localized soliton-like wave with the peak-value 
at m = 0. The width and peak-value of the wave decreases and increases respectively with 
increasing the value of \. This is obtained in all cases discussed here. 

III. RANDOM SYSTEMS 

The ballistic motion of an initially localized particle in a one-dimensional nonlinear ran- 
dom binary alloy has been observed recently. |i| Here we also show the ballistic motion of 
a particle in completely random nonlinear systems. The random nonlinear systems are 
characterized by random distribution of the nonlinearity parameter, \ with the values 
< Xm < Xmax- All the site-energies are assumed to be zero. The initial excitation is 
applied at the zeroth site. Furthermore, we assume that the value of Xo is Xmax- The MSD 
is defined as 

oo 

<m 2 >= rn 2 \c m (t)\ 2 (3) 

m=— oo 

with the initial condition |c m (0)| 2 = 5 m> o- After some initial transient behavior the speed 
(V < m 2 >/t) of the particle is setteled down to a constant which depends on Xmax- The 
speed of the particle decreases with increasing Xmax- Around the critical value of Xmax the 
speed decreases drastically and beyond this region we observe the slow decay of the speed. 
The untrapped portion of the probability leads to the ballistic motion of the particle above 
the critical value of Xmax- 13 To obtain the critical value of Xmax we have plotted the time 
averaged probability of the initial excitation site for different realizations in Fig. 8. The 
critical value of Xmax are found to be in the range of ~ 3 and ~ 3.5. Beyond the critical 
value, < Pq > in general increases very fast with increasing x but with a certain degree 
of sample to sample variation. Within the region of % ~ 3 and % ~ 4.5 we obtain a large 
deviation in < P > f° r different realizations. Beyond this region < P > increases slowly 
as Xmax goes up and the probability of all other sites decreases. In this limit the random 
nonlinear lattice to a good approximation can be replaced by a perfect lattice with a single 
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nonlinear defect of same strength (xmax)- This is also true for perfect nonlinear system. It 
can be shown that beyond the transition region of \ the speed of the particle in random 
nonlinear system, perfect nonlinear system and the single nonlinear impurity problem are 
almost equal. We study next the effect of nonlinearity in the superdiffusive motion of the 
random-dimer model (RDM). 

The RDM is a binary alloy containing two types of atoms with site energies e a and e&. 
The restriction on the randomness is that the site energy e a appears in a pair which is called 
a dimer. So, the system is the random distribution of the dimer and other site energy e&. If 
we assume e& = and V = 1 then we obtain ~ x/iV number of nonscattered states around 
p4| Here, N is the length of the sample. The MSD goes as < m 2 >~ t 3 ^ 2 . |H| This 



e n , 



is obtained only when |e Q | < 2. We now study the effect of nonlinearity on the transport 
properties of the RDM. We assume that all the sites have equal nonlinearity strength and it 
is x- It should be noted that by adding the nonlinearity the effective site energies are altered. 
Thus initially the dimer correlation is distorted. As the system contains escape probability, 
the distortion of the dimer correlation decreases with increasing time. So, we expect the 
superdiffusive behavior in the MSD. This is exactly obtained (see Fig. 9). The exponent 
of MSD is ~ 1.5. We also study the time averaged probability of the initial excitation site 
which is one of the dimer sites. We obtain a sharp self-trapping transition at different values 
of Xcr2 which increases as e a increases. Thus, by increasing the site energy we can alter the 
value of Xcr2- For negative values of e a the value of Xcr2 remains almost constant. If the 
initial excitation is applied at the site where dimer is absent we obtain the opposite behavior. 
That is, for positive values of e a the Xcr2 almost remain constant and for negative values of 
e a, Xcr2 changes significantly. Above Xcr2 we also obtain the superdiffusive behavior in MSD. 



IV. SUMMARY 

We have studied the dynamics of quasiparticles in different kinds of nonlinear systems. 
We first studied the dynamics of the quasiparticle in a cluster of nonlinear impurities em- 
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bedded in a perfect linear host lattice. The initial excitation is applied either at the end-site 
or at the middle-site of the nonlinear cluster. In both the cases we studied the time averaged 
site-probabilities and the particle propagation. In the former case we observed the cluster- 
trapping transition at Xcri due to the localization of the quasiparticle within the nonlinear 
cluster. For x > Xcri the amplitude-transition in the probabilities of the nonlinear sites 
is obtained showing that escape probability through linear sites is analogous to a damping 
term in oscillator equation of motion. The absence of any well defined transition in the 
amplitude of the site-probabilities beyond Xcr2 indicates the self-trapping at this value of x- 
This value of Xcri is ~ 4.23 for any size of the cluster of the nonlinear sites. This clearly 
indicates that a single mechanism is responsible for self-trapping. For relatively large size 
cluster we observed that the localization of the quasiparticle occurs in the cluster in the 
form of soliton-like wave extending over a few sites. This cluster-localization starts between 
X = 3.1 and 3.2. In the asymptotically large size of the nonlinear cluster we also find that 
soliton-like wave moves along the cluster but after some time it localizes in a region. This 
localization time and the trapping region depends on x- When the initial excitation is ap- 
plied at the middle site of the cluster, the cluster-localization also occurs. The critical value 
in general increases with increasing the size of the cluster. However, here we do not find any 
sharp self-trapping transition as well as the transition in the amplitude of site-probabilities. 

The MSD in random nonlinear system shows the ballistic motion and the speed decreases 
significantly in the transition regime. The initial excitation is applied in this case at the 
site with maximum nonlinearity. The time averaged probability of the initial excitation 
site is also studied for different realizations. In each cases we obtained rapid increase in 
< Po > within a range of x- It should be noted that beyond this region of x the dynamics of 
the quasiparticle in the single nonlinear impurity system, in random nonlinear systems and 
in the perfect nonlinear system are similar. The effect of nonlinearity on the superdiffusive 
behavior of the RDM has also been studied. The exponent of the MSD remains almost same. 
With increasing the strength of the nonlinearity the prefactor of the MSD decreases with 
the increase in x- The self-trapping transition also occurs in this case. The self-trapping 
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value of x increases with increasing e a provided \ an d e a both are positive quantity. For 
negative values of e a , this value of x remains almost constant with increasing x- This aspect 
will be studied latter. 
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FIGURES 

FIG. 1. Time averaged probability at the initial excitation site (< Po >) as a function of x for 
different size (n) of nonlinear cluster. The initial excitation is applied at the end-site (m = 0) of 
the cluster. 

FIG. 2. Time averaged probability of initially occupied (< P$ >) and unoccupied site (< Pi >) 
of the nonlinear dimer embedded in a host lattice as a function of x- The inset shows details of 
the < Pq > and < P\ > with x hi the fluctuation regime. The initial excitation is applied at the 
end-site (m = 0) of the dimer. 

FIG. 3. Plot of site-probabilities of the dimer sites embedded in a host lattice as a function of 
time (t) for different values of nonlinearity parameter x as follows: (a) x = 3-72; (b) x = 3.77; (c) 
X = 4 and (d) x = 4.4. In all these cases the initial excitation is applied at the end-site (m = 0) 
of the dimer. 

FIG. 4. Electronic probability propagation profile as a function of time (t) for x = 3.51. Here, 
n = 30 and the initial excitation is applied at the end-site (m = 0) of the nonlinear cluster. 

FIG. 5. Same as Fig. 4 but n — > oo and x = 3-6. 

FIG. 6. Time averaged probability of the initial excitation site (< Po >) as a function of x f° r 
different values of n. The initial excitation is applied at the middle-site (m = 0) of the nonlinear 
cluster. 

FIG. 7. Time averaged probability (< P m >) of the nonlinear sites of the cluster of size n = 5 
embedded in a host lattice as a function of %. The initial excitation is applied at the middle-site 
(m = 0) of the cluster. As the system is symmetric around the zeroth site < P_i >=< Pi > and 
< P-2 >=< P 2 >. 

FIG. 8. Time averaged probability of the initial excitation site (< Po >) as a function of Xmax 
for different realizations of random nonlinear system. 
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FIG. 9. Plot of < m 2 > /t 3//2 as a function of time (t) for the RDM with different values of %. 
Here, e a = 1. 
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